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Mass accretion onto black holes releases energy in the form of radiation and outflows. Although 
the radiative flux cannot substantially exceed the Eddington limit, at which the outgoing radiation 
pressure impedes the inflow of matter, it remains unclear whether the kinetic energy flux is 
bounded by this same limit. Here, we present the detection of a radio-optical structure, powered 
by outflows from a non-nuclear black hole. Its accretion disk properties indicate that this black 
hole is less than 100 solar masses. The optical-infrared line emission implies an average kinetic 
power of 3 x 10 40 erg second -1 , higher than the Eddington luminosity of the black hole. 

These results demonstrate kinetic power exceeding the Eddington limit over a sustained period, 
which implies greater ability to influence the evolution of the black hole's environment. 


F or an accreting black hole (BH) of mass M, 
the Eddington luminosity [L Edd ~ 1.3 x 
10 38 (M/M©) erg s -1 ; where M© is solar 
mass] is the maximum photon power achievable 
for spherically symmetric, radiatively efficient ac- 
cretion, and the Eddington accretion rate is the 
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corresponding mass accretion rate. However, BHs 
can also release energy through mechanical chan- 
nels (kinetic energy of jets and winds). The balance 
between radiative and kinetic output depends 
partly on the mass accretion rate. For accretion 
rates -0.03 to 0.5 times the Eddington rate, the 
BH power output is entirely dominated by ther- 
mal disk emission without a jet (1-3), whereas 
for rates approaching or exceeding the Eddington 
rate, we expect radiatively driven outflows, pos- 
sibly coexisting with faster, collimated jets (4, 5). 

There is circumstantial evidence that non- 
isotropic photon luminosity can mildly exceed 
the Eddington limit (3, 6-8), but can BHs with 
supercritical accretion have powerful jets, and 
can the mechanical power Pj exceed L E dd for an 
extended period of time? Empirical evidence of 
super-Eddington jet power was presented for a 
sample of radio-loud quasars (9, 10). However, 
ambiguities in the BH masses and in the con- 
version from 151-MHz radio flux to jet power 


(11-15) make these results uncertain. In the local 
universe, shock-ionized gas around ultraluminous 
x-ray sources (ULX bubbles) provides strong evi- 
dence of sustained and powerful mechanical 
output from sources that are thought to have high 
accretion rates (16-22). Typical ULX bubbles 
have sizes of -50 to 300 pc and require a me- 
chanical input power of -10 39 to 10 4 ° erg s -1 for 
-10 5 years. However, there are no BH mass mea- 
surements for such ULXs: If they are powered 
by -10 M© BHs, their radiative and/or mechan- 
ical power would, in some cases, exceed the 
Eddington value by a factor of 10, but if they 
are powered by BHs of mass >100 M© (known 
as intermediate-mass BHs), they may not exceed 
the Eddington limit. In the Milky Way, the most 
energetic BH bubble is probably SS 433/W50, with 
an average mechanical power of -10 39 erg s -1 over 
an age of a few 10 years (23-28). 

Here, we report the detection of a powerful 
microquasar in the inner disk of the spiral gal- 
axy M83 (NGC 5236), located at a distance of 
4.61 ± 0.26 Mpc (29). We are conducting a multi- 
wavelength study of x-ray source populations in 
M83. Our observations include 5.5- and 9.0-GHz 
images taken with the Australia Telescope Com- 
pact Array (ATCA) on 28 to 30 April 201 1, for a 
total observation time of 36 hours. We also used 
archival data of M83 taken with the Australian 
Long Baseline Array (LB A) on 13 May 2005, for 
a total of 9 hours, at 2.3 GHz (30). In the infrared 
(IR)-optical bands, we used archival Hubble Space 
Telescope/Wide Field Camera 3 (HST/WFC3) 
images of the inner disk of the galaxy in several 
broad- and narrow-band filters (data taken on 
20 to 26 August 2009). We observed M83 with 
the Chandra X-ray Observatory/ACIS-S camera 
10 times between December 2010 and December 
2011, as part of a large program to study x-ray 
source populations of a star-forming galaxy; the 
total exposure time was 729 ks. In addition, we 
used an archival 50-ks x-ray observation of M83 
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Fig. 1. Location of MQ1 in the ^ 0 

inner region of the M83 galaxy. 

(A) HST/WFC3 image of the central 
region of M83 in the F657N filter, 
with superposed 9.0-GHz surface 
brightness contours, from our ATCA 
data. Contour levels are 0.10, 0.18, 

0.40, 0.78, 1.32, and 2.0 mjy per 
beam. The beam size is 2".3 x 0".9 
at a position angle of 173°. The 
root mean square noise level in 
source-free regions of the map is 
0.025 mjy per beam, and the peak 
intensity of MQ1 is 0.9 mjy per 
beam, but this includes a small con- 
tribution from the surrounding ex- 
tended emission. Labels mark the 10 '' _ 220pc w i 

location of the microquasar MQ1 4 0 

and of a nearby SNR used for astro- 

metric alignment. The blue circle (radius 0".6) marks the nucleus of M83. (B) Chandra/ACIS x-ray color image of the same field (red, 0.35 to 1.1 keV; green, 
1.1 to 2.6 keV; blue, 2.6 to 8.0 keV). 
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Fig. 2. X-ray hardness and spectral properties. (A) X-ray hardness-versus- 
intensity diagram of MQ1 from our sequence of Chandra observations. The 
number next to each data point indicates the order in which the observations 
were taken (table SI). The blue line denotes the expected location of a disk- 
blackbody spectrum with temperature /r B 71 n varying (from top to bottom) from 
0.80 to 0.50 keV in steps of 0.05 keV. Dashed and dash-dotted lines mark the 



expected hardness of sources with the same absorbing column and power-law photon 
indices r = 2.5 and 3, respectively: The observed values are clearly inconsistent with 
the low/hard state (r < 2) at all epochs, ct s _1 , counts per second. (B) Combined x-ray 
spectrum from the Chandra ObsID 12994 (6th observation) and 14342 (11th 
observation), when the source was brightest. The spectrum is well fitted with an 
absorbed disk-blackbody model, with parameters typical of stellar-mass BHs (table S3). 



Fig. 3. Optical-IR line emission. (A to D) HST/WFC3 image of the microquasar in four continuum- 
subtracted bands. Clockwise from top left: [Fell] X = 1.644 pm; [Sll] M6716, 6731; [Olll] X5007; and 
Ha + [Nil] AA6548, 6584. The blue circle indicates the centroid of the ATCA/Chandra source and has 
a radius of 0".2. 


obtained with Chandra/ACIS-S in April 2000 
(31, 32). 

Our ATCA images show a luminous radio 
source at right ascension (RA) =13 hours 37 min 


1.123 s (± 0.005 s), Dec. =-29°51'52".17 (± 0".10) 
(J2000), ~5" northeast of the optical nucleus of 
M83 and just outside the intensely star-forming 
region that surrounds it (Fig. 1 and fig. SI). It has 


a flux density S 5 _ 5 = 1.4 ± 0.1 millijansky (mJy) 
and S 9 q = 0.80 ± 0.05 mJy at a frequency (v) of 
5.5 and 9.0 GHz, respectively, which corresponds 
to a radio luminosity ~2 x 10 35 erg s -1 , similar to 
those of the two most powerful microquasar 
bubbles known to date, IC 342 X-l and NGC 
7793-S26 (19, 20, 22). At 5 GHz, this radio 
source is ~30 times more radio-luminous than 
SS 433/W50 (24, 26, 27). Henceforth, we identify 
this source as MQ1. Its steep spectral index a ~ 
1.1 (S v ~ O implies optically thin synchrotron 
emission, possibly from an aging population of 
electrons. MQ1 is unresolved in the ATCA maps, 
implying a diameter < 1" ~ 22 pc. In the LBA 
map, an unresolved source (size < 2 pc) is detected 
atRA= 13 hours 37 min 1.136 s (± 0.003 s), 
Dec. = -29°51'52".15 (± 0".02) (J2000) with a 
flux density of S 2 .3 = 0.22 ± 0.04 mJy at 2.3 GHz, 
a factor of ~10 to 15 times less than the total 
expected flux density at that frequency, based on 
the ATCA data. This indicates that -90% of the 
emission is extended over a region of at least a 
few parsecs in diameter (much larger than the 
LBA beam but smaller than the ATCA beam). 
The LBA source is offset to the east of the ATCA 
centroid by ~0".2 ± 0".l (fig. S2). 

We interpret the radio source as emission from 
a radio bubble or radio lobes (that is, from an ex- 
panding structure of shocked plasma) inflated by 
the mechanical power of the BH; synchrotron- 
emitting electrons are produced where the jet 
heads hit the interstellar medium (ISM). Similar 
structures have been found in other non-nuclear 
BHs: notably, the ring around Cyg X-l (33), the 
“ears” of SS 433 (34), and the bubbles around 
several ULXs (35). We speculate that the LBA 
radio source arises from the eastern hot spot, 
where the head of a jet strikes the ISM. The 
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nondetection of a compact core in the LBA 
(flux density S 2 . 3 <0.1 mJy per beam) does not 
imply that the jet is currently turned off, as we 
estimate [scaling from the power and luminosity 
of galactic microquasars, e.g., (36)] that the ra- 
dio flux density of the MQ1 core would always 
be less than or equivalent to a few microjansky, 
even for an efficient jet production at high ac- 
cretion rates. 

Our Chandra study shows a luminous point- 
like source, coincident with MQ1 to within the 
x-ray positional uncertainty of 0".3, whose char- 
acteristics provide an estimate of the BH mass. 
The x-ray flux and hardness vary from one obser- 
vation to another (Fig. 2 and fig. S3) and also 
within observations (fig. S4). This variability mles 
out emission from shocked gas; the source must 
be an accreting compact object. The x-ray spec- 
frum is always consistent with thermal disk emission 
(2); it is never consistent with the power-law- 
dominated low/hard state. When the source is 
brightest [Chandra X-ray Observatory identifica- 
tion number (ObsID) 12994 and ObsID 14342], a 
standard disk-blackbody model fit [tbabs x 
tbabs x diskbb in XSPEC version 12.6.0 (37)] 
gives a characteristic peak temperature k B T m ~ 
0.7 keV (where k B is the Boltzmann constant), 
apparent inner-disk radius r in (cos0) 1/2 ~ 50 km 
(where 0 is the viewing angle: face-on for 0 = 0), 
and unabsorbed 0.5- to 8-keV luminosity ~ 7.4 x 
10 37 (cosO) 1 erg s -1 . The inferred BH mass de- 
pends on the viewing angle and the BH spin: M~ 
(\.\9c 2 IG)(\la)r m (here, c is the speed of light 
and G is the gravitational constant) (38), where 
a = (c 2 /G)(Ri sco /M) = 6 for a nonrotating BH and 
1.24 < a < 6 for astrophysically plausible 
spinning BHs (39, 40), and R isco is the innermost 
stable circular orbit. From our standard disk- 
blackbody fit, we obtain M ~ (42.6j 9 12-8 )(l/a) 
(cosO) 1 / 2 M q . Thus, for arbitrary midrange val- 
ues (cosO ~ 0.5, Kerr parameter a ~ 0.5, implying 
a ~ 4), the spectral model, fitted temperature, 
inner-disk radius, and luminosity are similar to those 
found in galactic stellar-mass BHs (M ~ 10 M©) 
accreting at -10% of the Eddington rate, which 
supports the stellar-mass BH identification for 
MQ1 (Occam’s razor). For the most extreme (less 
likely) choice of parameters in their physically 
acceptable range (0 = 85°, a - 0.998), the best- 
fitting BH mass estimate is M~ 115 M©. 

There is an optical emission line source at the 
position of MQ1 that was previously identified 
with HST (41) and classified as a supernova rem- 
nant (SNR) within the nuclear region. Closer 
inspection of the optical source in the continuum- 
subtracted Ha + [Nil] M6548, 6584, [Fell] 
XL 644 pm, and [SII] XX6716, 6731 A images 
(Fig. 3) shows that it is resolved into a stronger 
central source with two knots or lobes at either 
side (a structure reminiscent of SS 433/W50). 
The sum of the projected distances between the 
side knots and the central peak is ~0".6 ~ 13 pc. 
There are no other bright stars or sources of op- 
tical continuum emission at this location (fig. S5); 
thus, the extended optical source is a physical 


structure rather than a random alignment of three 
or more unrelated emission line objects. We used 
nearby sources detected in more than one band 
to refine the relative astrometry of the HST image 
with respect to the radio coordinates, and we con- 
clude that the optical core, the ATCA centroid, and 
the x-ray position are coincident within ~0". 15. The 
eastern knot of Ha emission may be associated 
with the compact LBA source (fig. S2). We sug- 
gest that the optical structure represents the core 
and hot spots and/or lobes of the microquasar. 

Optical line emission from the shocked gas is 
a reliable proxy for mechanical power. For a fully 
radiative shock expanding into the ISM with ve- 
locity v s , the total radiative flux equals the flux of 
mechanical eneigy through the shock (35, 42). Stan- 
dard bubble theory (43) shows that the bolometric 
photon luminosity L T is related to the long-term- 
averaged mechanical power as L T = (21/ll)P h 
with the rest of the power going into the kinetic 
energy of the expanding shell of swept-up gas and 
into the thermal energy of the shocked gas between 
the reverse shock and the swept-up shell. The frac- 
tion of bolometric radiative luminosity emitted in 
a given line (including both the shock and the 
precursor contributions) can be computed from 
shock- and photoionization codes such as Model- 
ling And Prediction in Photo-Ionized Nebulae and 
Gasdynamic Shocks III (MAPPINGS III) (44). 

Hp or Ha are the lines most often used to 
estimate microquasar jet power (19). However, 
in this case both lines (especially Hp) are affected 
by high optical extinction; the Ha flux from HST 
imagery includes an uncertain but substantial con- 
tribution from [Nil] XX6548, 6584; and part of 
the Balmer emission may come from x-ray photo- 
ionization due to the central source. Thus, we used 
instead the flux in the IR line [Fell] X = 1.644 pm 
as a proxy for the mechanical power. The [Fell] 
line is a powerful diagnostic tool for shock-heated 
gas in high-density regions (45-48). 

We measured the [Fell] line flux (table S2) 
from the HST/WFC3 IR-camera image in the 
F164N filter after subtracting the continuum 
emission using the F160W image. The apparent 
(not corrected for extinction) [Fell] luminosity is 
~1.1 x 10 37 erg s- 1 , one order of magnitude greater 
than that observed from SNRs in the spiral arms. 
We adopted an extinction A v ~ 3.9 magnitude 
(mag), as inferred from the x-ray column density 
[n H = (8.6^5) x 10 21 cmT 2 ] using the conver- 
sion relation of (49), with a likely uncertainty 
range A v - 3.8 to 5 mag. Taking into account 
the corresponding IR extinction at 1.64 pm 
(A 1 ,64|j.m) ~ 0.70 mag, the emitted luminosity 
Z[Feii] = (2.1 + 0.2) x 10 37 erg s _1 (table S2). 

To estimate the mechanical power, we assume 
a preshock atomic number density n 0 -100 cm -3 , 
typical of the inner disk-nuclear region of M83 
(41, 50, 51). Then, from MAPPINGS m models 
(44), and assuming equipartition between mag- 
netic and gas pressure, we derive that T [FeII ] = 
(0.7 ± 0.1) x 10 -3 Pj over a wide range of plau- 
sible shock velocities (fig. S6). We conclude that 
Pj ~ 3 x 10 40 erg s -1 > L Edd (from the BH mass 


constraints). We can derive other physical param- 
eters using a well-known self-similar expansion 
solution for jet-driven bubbles and lobes (43, 52). 
The characteristic radius R (distance between core 
and hot spots) scales as R - 0.76Pj 1/5 t 3/5 p 0 ~ 1/5 
where p 0 - p m p n 0 (m p is the proton mass and 
p the mean atomic weight). Taking Pj ~ 3 x 
10 40 erg s- 1 , R - 0".3 ~ 6.7 pc (neglecting projec- 
tion effects), and p = 1.38, we infer a character- 
istic expansion age t - 4.9 x 10 11 s - 16,000 years. 
The expansion velocity v s = dR/dt = (3/5 )R/t - 
250 km s' 1 . Taking n 0 - 50 to 200 cm -3 does 
not substantially change the power but gives an 
age range -13,000 to 20,000 years and a veloc- 
ity range v s ~ 200 to 310 km s -1 . Mechanical 
power and shock velocity are similar to those 
measured in the powerful microquasar NGC 
7793 S26 (19, 20, 53), but the age of MQ1 is a 
factor of 10 younger and the ISM density two 
orders of magnitude higher. 

Both the optical-IR and the radio luminosity 
come from the bubble and/or lobes; therefore, 
they trace the time-averaged mechanical power. 
The x-ray luminosity traces the current accretion 
rate. It is likely that during the Chandra obser- 
vations, the system had a sub-Eddington luminos- 
ity and a weak jet or none at all. It is also likely 
that in the past, the system had super-Eddington 
luminosity and jet power: It would have looked 
like ULXs such as Holmberg IX X-l or NGC 
1313 X-2 (luminosity L x ~ 2 x 10 40 erg s -1 ), sur- 
rounded by ionized bubbles with mechanical power 
P] - 10 40 erg s' 1 (35). If the x-ray core switches off 
or fades, the surrounding radio-optical bubble re- 
mains bright. 

For MQ1, we have been able to constrain the 
BH mass from x-ray spectral analysis because 
the BH is currently in a disk-dominated state. 
None of the other microquasars with similar me- 
chanical power has a reliable constraint on its 
BH mass, because none of them has been ob- 
served in this state. Combined with an accurate 
determination of time-averaged mechanical power 
from shock-ionized line luminosities, this enables 
a robust estimation of the (mechanical) Eddington 
ratio. Thus, we have identified a system that has 
had (on average) a super-Eddington mechanical 
power for -2 x 10 4 years (most likely, a few 
times above Eddington). MQ1 provides empirical 
evidence that BHs can sustain super-Eddington 
mechanical power for an extended period of time. 
During those phases, BHs inject more energy into 
the surrounding medium than would be inferred 
on the basis of their Eddington limits. The exis- 
tence of super-Eddington mechanical power is 
important for the modelling of jet production 
mechanisms, as well as the evolution of fast- 
growing supermassive BHs in the early universe 
and their effect on their host galaxies. 
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Large-Amplitude Spin Dynamics 
Driven by a THz Pulse in Resonance 
with an Electromagnon 

T. Kubacka, 1 * ]. A. Johnson , * 2 M. C. Hoffmann , 3 C. Vicario , 4 S. de Jong , 3 P. Beaud , 2 S. Griibel , 2 
S.-W. Huang , 2 L. Huber , 1 L. Patthey , 4 Y.-D. Chuang , 5 J. J. Turner , 3 G. L. Dakovski , 3 W.-S. Lee , 3 
M. P. AAinitti , 3 W. Schlotter , 3 R. G. Moore , 6 7 C. P. Hauri , 4,7 S. M. Koohpayeh , 8 V. Scagnoli , 2 
G. Ingold , 2 S. L. Johnson/t U. Staub 2 t 

Multiferroics have attracted strong interest for potential applications where electric fields control 
magnetic order. The ultimate speed of control via magnetoelectric coupling, however, remains 
largely unexplored. Here, we report an experiment in which we drove spin dynamics in multiferroic 
TbMn0 3 with an intense few-cycle terahertz (THz) light pulse tuned to resonance with an 
electromagnon, an electric-dipole active spin excitation. We observed the resulting spin motion 
using time-resolved resonant soft x-ray diffraction. Our results show that it is possible to 
directly manipulate atomic-scale magnetic structures with the electric field of light on a 
sub-picosecond time scale. 


D ata storage devices based on ferromag- 
netic or ferroelectric materials depend 
strongly on domain reorientation, a pro- 
cess that typically occurs over time scales of 
several nanoseconds. Faster reorientation dy- 
namics may be achievable by using intense 


electromagnetic (EM) pulses (7). The EM pulses 
can couple to magnetism either indirectly via 
electronic excitations (2) or directly via the 
Zeeman torque induced by the magnetic field 
(3-5). Direct excitation has the advantage of 
minimal excess heat deposition but requires 


frequencies in the 10 10 to 10 12 Hz range. The low 
magnetic field strength of currently realizable 
THz-frequency EM sources poses a formidable 
challenge for such schemes. 

Thanks to the coexistence of different ferroic 
orders, multiferroics offer new routes to domain 
control ( 6 ). Particularly strong coupling between 
the ferroelectric and magnetic order exists in 
single-phase frustrated magnets, where noncol- 
linear spin structure drives ferroelectricity as a 
result of weak relativistic interactions (7-9). 
Consequently, the magnetic order can be con- 
trolled by application of an electric field (10-13). 
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